Behavior of the two-dimensional Ising model 
at the boundary of a half-infinite cylinder 



Yvan Saint- Aubin* Louis-Pierre Arguin^and Hassan Aurag-^ 

September 22, 2009 

Abstract 

The two-dimensional Ising model is studied at the boundary of a half-infinite cylin- 
der. The three regular lattices (square, triangular and hexagonal) and the three regimes 
(sub-, super- and critical) are discussed. The probability of having precisely 2n spin- 
flips at the boundary is computed as a function of the positions fcj's, i = 1, . . . ,2n, 
of the spinflips. The limit when the mesh goes to zero is obtained. For the square 
lattice, the probability of having 2n spinflips, independently of their position, is also 
computed. As a byproduct we recover a result of De Coninck showing that the lim- 
iting distribution of the number of spinflips is Gaussian. The results are obtained as 
consequences of Onsager's solution and are rigorous. 



1 Introduction 

The present paper gives a rigorous description of the spin configurations seen at the end of a 
half-infinite cylinder covered by the Ising model on square, triangular or hexagonal lattices. 
Both the discrete case and its continuous limit are made explicit. Section 2, devoted to 
the square lattice, describes the boundary behavior by giving, for a fixed number 2n of 
spinflips, the probability distribution as a function of the positions of the spinflips. The 
three regimes, subcritical, critical and supercritical, are discussed. It should be stressed that 
these probabilities are not correlation functions, even though they can be used to calculate 
them. In Section 3 we compute, for the square lattice, the distribution of the random variable 
(number of spinflips) / (number of sites at the boundary) at the critical temperature. Section 
4 extends the results of Section 1 to the triangular and hexagonal lattices. The results (and 
the methods to obtain them) are simple, though non-trivial, and the thermodynamical limits 
are shown to depend only on the behavior of smooth functions defined on [— vr, vr] for the 
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square lattice and on [— 7r/2,7r/2] for the other two. Only their behavior at zero and at the 
extremities on these intervals play a role. 

The calculation presented here could have been done by several methods and for other 
orientations of the lattices. We chose the most classical technique, the one based on the 
transfer matrix; in the present case, this method turns out to be simple and scientists with 
expertise in neighboring fields will be able to follow the argument with minimal investment. 
It is impossible to make a short overview of more than sixty years of developments on the 
Ising model. There are two results however that we emphasize since they are directly related 
to ours. Abraham [1] proved that the limit distribution, with a non-trivial scaling, for the 
magnetization of the Ising model on a square lattice is Gaussian at the edge of a half-infinite 
cylinder and at criticality. Along the same line, De Coninck [6] showed that in the limit 
the joint distribution for the magnetization and the energy (which is a linear function of the 
number of spinflips) is also Gaussian. The two proofs use the transfer matrix method to 
directly compute the characteristic function of the variables. We recover here De Coninck's 
result for the number of spinflips using the distribution in the discrete case computed in 
Section 2 and a combinatorial lemma found in the Appendix. We obtain explicit expressions 
for both the mean and the variance of the limiting distribution. 

Finding the relative weights of boundary configurations as a function of the number and 
locations of spinflips might seem only a mildly interesting exercise. A word of explanation is 
therefore in order. With the invention of the stochastic Loewner equation (SLE) by Schramm 
[19] and proofs by Smirnov [21] that percolation and the Ising model are conformally invariant 
in the limit when the mesh goes to zero, new rigorous tools have been available to probe 
critical phenomena, and with these tools, new observables have been introduced. Suppose 
that boundary conditions are imposed at the extremity of the half-infinite cylinder as follows. 
Let < < < ^3 < ^4 < 277 be four angles. They define four intervals along the boundary 
and we suppose that the sign of Ising spins are constant on each, but alternate from one 
interval to the next. Such boundary conditions force interfaces, that is contours between 
constant-sign clusters, to intersect the boundary at the 6''s, and only there. Since these 
interfaces cannot cross, the interface starting at 6i must end at 62 or 6'4. One can therefore 
ask what is the probability that the interface starting at 61 goes to 62- At least five groups in 
the last eight years computed this new "observable" , and they all had to solve the problem 
that is the subject of the present paper. The first two groups, two of the present authors 
[2], and Bauer, Bernard and Kytola [3], used conformal field theory (CFT) to do it, even 
though the latter group was actually interested in devising a way to define multiple SLE 
processes. The three other groups used purely SLE methods to obtain the result [201 S IS] . 
We were not able to compute directly from the lattice models the probability just described. 
This is why we turned in [2] to CFT to do it. And this is why the present paper limits 
itself to computing the relative weights as functions of the number and locations of spinflips, 
disregarding how the interfaces join the spinflips. Despite its limitations, the calculation 
has several redeeming features. First, it is done from first principles and relies on classical 
methods. Second, contrarily to the results obtained from SLE or CFT, it gives an explicit 
result for any number of spinflips, not only for four. Third, it allows for subcritical and 
supercritical regimes to be studied. Finally, it provides information on a mesoscopic scale, 
namely on the distribution of the number of spinflips at the boundary. We note also that 
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the continuous distribution for the positions of four spinflips at the boundary (cf. equation 
(IT^ ) was an important element in the computation just described [2]. 

2 Behavior at the boundary for the square lattice 
2.1 Notation 

Let a : {1, 2, ... , m} {+1, —1} be a configuration along the circular extremity of a half- 
infinite cylinder covered by the square lattice. The number m of sites on this circle is taken 
to be even. It is convenient to encapsulate the information contained in a in the following 
data: the value s G {+1, —1} of a at 1 and the positions ki G {1, 2, . . . , m}, 1 < z < 2n, 
of spinflips. A spinfiip in a occurs at ki if a^-^i = —cr^- with ctq = and o"i = (Xm+i by 
definition. We choose 1 < ki < k2 < ... < k2n < fn. The number n can also be seen as the 
number of maximally connected stretches of +-spins along the boundary. We identify the 
functions a with their data (s; ki, ^2, . . . , k2n). 

If cr and a' are two configurations on contiguous circles along the cylinder, the trans- 
fer matrix T : C^" C^™ is given by its matrix elements T^j.^ = exp(i^ YlT=i '^fc'^fc+i + 
^ XlfcLi '^feCTfe)- The constant u is the product of the coupling constant of the Ising model 
taken here to be isotropic, with the inverse temperature measured in units that make the 
Boltzmann constant unity. We shall refer to u as the inverse temperature. The critical 
temperature is defined by sinh2z/c = 1. When the Ising model is on a torus, there is some 
freedom in the choice of the transfer matrix. Here this choice is unique as the transfer ma- 
trix must include the Boltzmann weights attached to the bonds between the m sites at the 
boundary and the bonds that tie them to the first inner circle. 

Note that all matrix elements are positive. By Perron-Frobenius theorem, the transfer 
matrix has real eigenvalues. Moreover its largest eigenvalue is non-degenerate. Let a; be a 
non-zero eigenvector corresponding to this eigenvalue and let c'^{ki,k2, ■ ■ ■ , k2n) be its com- 
ponents in the state basis. It is known that these components can be chosen such that they 
are all positive and that their sum is 1. With this choice, the probability of a given state a 

Pr(cr = (s; fci, /c2, • • • , k2n)) = c\ki, k2, . . . , /i;2n)/« (1) 

where 

m/2 

5Z 5Z c'{ki,k2,...,k2n). 

se{+l,-l} n—O I<ki<k2<...<k2n<m 

Note that the c^(ki,k2, ■ ■ ■ , A;2n)'s do not depend on s and this super- index will be deleted 
when it is appropriate. 

We use Thompson's notation [TS] in the rewriting of T in terms of tensorial blocks. The 
Pauli matrices = ( ? J ) , = ( ° "^^ ) , = ( o ~i ) , and 1 = ( q ? ) are used to define 
2m ^ 2''"-matrices = 1 ® ■ ■ ■ ® r* (g) • ■ ■ (g) 1, i G {1,2,3} and 1 < k < m, where the 
only non-trivial factor is at position k. The following operators (C^™ —>■ C^™) will also be 
used: pk = 'tI'tI ■ ■ ■ '^l-i'^l ^ ^"^^ '^k = '^I'^'l ■ ■ ■ '^l-i'^l where again k G {1,2, ...,m}. The 
linear combinations = ^{pk + iT^k) and al = \{pk — i^^k), k G {l,2,...,m}, satisfy 
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o-fcCffc' + a|,,afc = and a^ak' + Qfc'Qfc = 02™. Finally their Fourier coefficients are given 

by 

Vg = ^^ ^'''''"'k and ^J = ^ I] ^'"'4 

^ l<k<m ^ l<k<m 

for q G Qm, the set of phases of the m-roots of —1: Qm = ~y + 1 < j < yj- 

Using this notation, the transfer matrix is 

T = (2sinh2z/)^ exp(2z/7riPmP) exp I -iu ^ vr^+ipfc | exp I iz/* ^ n^Pk 

\ l<k<m-l / \ l<k<m 

where P = Y[i<k<m'^k^ ~ -'-a'"' is the operator ffipping a configuration a = {s^ki) into 
{—s]ki) and u* is defined implicitly by sinh 2z/ sinh 2i^* = 1. Since the Boltzmann weights 
are invariant under the flip of all spins, P commutes with the transfer matrix and T and P 
can be diagonalized simultaneously. The eigenvector u belongs to the +-eigensubspace of P. 

The eigenvectors of T are constructed with creation and annihilation operators ^| and 
^q,q G Qm, obtained from the rjg and by orthogonal transformations. They are = 
rjq cos (pq + rj'ig sin 0^ and = rj^g cos (pg — 77^ sin 0g with 

(tanh 2z/ + sinh 2z/) sin q 



tan0q 



1 — 2 sinh 2z/ cos q + a/ (cosh2z/ — tanh2z/(cosg + l))(cosh2z/ — tanh2z/(cosg — 1)) 



Though our choice of the transfer matrix, and therefore our expression for tan (pg, are different 
from those used in [22], the reader will find in that reference the method to obtain this ex- 
pression. The eigenspace spanned by u is characterized algebraically as the one-dimensional 
kernel of the m operators C,q'- 

iguo = 0, g G Qm- (2) 

The other eigenstates of T in the maximal subspace V+ G V where P\v+ = 1 are obtained 
by acting with an even number of on the vacuum u). Because the C,q and anticommute 

like the pairs ak, a\/s and rjq, r]^q/s, a vector ^l^^^l^ ■ ■ ■ ■Cga,'^ is a non-zero eigenstate only if the 
gi's are distinct. Equation ([2]) is the one to be solved. 



2.2 The discrete and continuous cases n = 1 

Suppose a describes a configuration with a single stretch of — spins and a single stretch 
of +-spins. Then a = {s;ki,k2) and n = 1. This paragraph is devoted to computing 
Pr((T = (s; ki, k2)) in the discrete case, cf. equation ([3]), and in the limit, cf. Proposition [3l 
By translation invariance, it is sufficient to compute c^(l,fc). The operators flip the 
spin at position k and the operators pk and i^k+i flip all the spins from 1 to inclusively. The 
operators ^q are linear combinations of these and the component along cr-f (the configuration 
with only +-spins) in E,qOJ can therefore originate only from the action of S,q on a-f or on a 
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a of the form (— ; 1, /c). If we denote by {u)^ the component along v in the vector u (in the 
basis given by the cr's), then 

jp-i-K/A ( ^ 



and 



2<fe<m 



^ ^ 2<k<m 



Using C| = Cj and e*^*" = —1 for q e Qm, the vanishing of {^qUj)cr^ gives i Yli2<k<m ^ ^'''^ (-'-' ^) 
c-f cot(0q + 2)- Using the discrete inverse Fourier transform and translation symmetry, one 
gets: 

c^(A;i, k^) = E e^^^'^-^^) cot(0, + |). (3) 



q&Qn 



That c^(A;i, ^2) is a real number follows from c^{l, k) = c^(l,m + 2 — /c) and the fact that 
(pq is an odd function of q. The expression for cot(0q + |) is surprisingly similar to that of 



tan (pq-. 



d{v, q) = cot(05 + I) 

(1 — tanh 2y) sin q 



sinh 2v — cos q + •>/ (cosh 2i/ — tanh 2i/(cos + 1)) (cosh 2v — tanh 2i/(cos g — 1)) 

(4) 

We note as a curiosity the following simple relation: cot(0q + f) ~ tan 0g at = i^c- 
We gather in a lemma the elementary properties of diy^ q) that are needed for the limit. 

Lemma 1 The function d : M+ x M ^ R defined by d{i', q) — coi{(j)q + q/2) given above has 
the following properties: 

(i) for any v e (0, 00) , d has a simple zero at q — mr, n odd; 

(ii) around q — Q (or any q = 2n7r,n G Z^, the function has the following behavior 

ly > Vc (subcritical) d{v, q — Q) — Q and is smooth around this point 
ly — Vc (critical) d^u, q) has a jump at q = with 

\f2-\ = limg^o+ '^('^c, (l) = - linig^o- c^(z^c, 
V <Vc (supercritical) d{y^ q) has a simple pole at q = 0; 

(iii) outside the behavior along the lines q — 2mr stated in (ii), d{i', q) is analytic in q for 
all inverse temperatures v. 
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-1 + V2 




Figure 1: The three regimes of the function cot(0q + g/2): subcritical (dashed curve), critical 
(plain curve) and supercritical (dotted curve). 

Figure 12.21 draws a graph of these three regimes. 

Proof: Property (i) is obtained by direct evaluation. For (ii) and (in) note first that the 
argument of the square root appearing in the denominator vanishes only when = 
cosg ± 1. Because the Ihs has a single minimum at = z/c = ^arcsinhl and is equal to 2, 
the argument of the square root is zero for v = and q G 27rZ only. 

For V > Uc, the square root is positive since sinh2z/ > 1 > cosg, and the denominator 
never vanishes. The analyticity of the subcritical case follows. For < i>c, the extrema of 
the argument of the square root can be found to be at g G ttZ with minima at g = 2nn. The 
denominator at g = nvr, n odd, is equal to 1 + sinh 2z/ + a/ cosh 2z/(cosh 2^ + 2 tanh 2z/) and is 
therefore positive. To check whether it vanishes it is sufficient to evaluate it at its minimum 
g = 0. At this point the denominator is sinh2z/— 1 + a/(1 — sinh 2//)^. Since sinh2z/c = 1, we 
have that the denominator vanishes only when g = and ly < I'c- For u < Uc, the expansion 
around g = of the denominator is ^(cosh^ 2z/(l — sinh 2i/))~^ + O(g^). Because of the sin g 
in the numerator, the supercritical behavior described in (ii) and (Hi) is established. For 
u = z/c, the function d takes the simple form 

(V2-l)sing 

rf(z/c, g) - 



a/2(1 — cos g) + a/(1 — cosg)(3 — cosg) 

The leading behavior of its denominator around g = is \/ q'^ + 0{q^) + ^ + Oi^q"^) which 
accounts for the jump. □ 
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Before turning to the continuum result, we need a technical lemma that rewrites the limit 
of equation ([3D when d is smooth. Equation ^ shows that c^{ki, ^2) depends only on the 
ratio (^2 ~ ki)/m since iq{k2 — fci) = i{2j — l)7r(A;2 — ki)/m. The meanfingful thing to do 
is to obtain the limit of that expression when ^ = G Q fl (0, 1) is held fixed. For this 
fixed ratio, the m's involved in the limit will have to verify two conditions: (1) If ^ G Q 
is written as - for two relatively prime integers, then n must divide the m's used in the 
sequence. This ensures that ^ is actually of the form — ki)/m. The second condition is 
more technical and its meaning will be clear in Section SI In the case of the square lattice, 
the momenta q can be restricted to the interval [— 7r,7r]. For the triangular and hexagonal 
lattices, they lie in [— f , f ]• We introduce 7 G N such that the momentum range is [— ^5 
The second condition on m's is: (2) r := ^ mod n is constant. 

Lemma 2 (i) Let 2; : [0, ^] — M, 7 G N 6e such that z,z' and z" exist and are continuous 
on [0, Let ^ = f with gcd(p, n) = 1. Then 

Jim 2^si„ ((2, - ((2, - 1)^) = -— (.(0) - .ef)cos(rO)) (5) 



m/7 



lim 2 Vcos (i2j -l)^-)z (i2j - 1)^) = ^^.(^) (6) 
m^oo ^ ^2/ ^mJ sin^/2 



where the limits are taken over integer m's such that n and 7 divide m and such that r := m/'j 
mod n is constant. 

(ii) If z(^) = 0, the above results hold without the requirement that r be held constant. 

Proof: (i) We divide the range of j G {1, 2, . . . , into subsets of n consecutive elements: 
{1,2,..., n}, {n + 1, n + 2, . . . , 2n} and so on. There will be r elements left. Omitting these 
r terms for the present, we consider the sum 

r^l_i r "1 1 1 

'-n'y' n — 1 '-nj' n—1 

^^'■= H 5^e^(2(^'+'")+i)^/2^((2(j + /n) + l)7r/m) = ^ ^ e^(2,+iW2^(^^ ^ 2j7r/m) 
1=0 j=0 1=0 j=0 

(7) 



where qi = {2ln + l)7r/m and [x] stands for the integer part of x. By Taylor expansion 

^eij with \eij\ < s := supgg[o,2V7] 



z{qi + 2j7i/m) = z{qi) + "^z^qi) + ^^e^j with < s := sup groo^/^i \z"{q)\ < 00. The 



two sums X]"=o e*^^-''*'^''^^^ and X]"=o J'^*''^'''''^^^^^ easily calculated. At = |, the first 
vanishes and the second is „. ^am - The sum (ITl) is then 



[— — 1 

I- J n—1 



/=o i=o ^ 



[—1-1 [—]-!», 1 

n-K 1 X ^ , 27r \ ^ X ^ 

-■ 1^ z\qi) + — 2^ 

imsmU 2 ^ — ^ m"^ ^ — ^ ^ — ^ 

' i=o «=o j=o 
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2gi{2j+l)e/2g . 



The second sum goes to zero since its absolute value is smaller than l-ji^v^sj^m. The first 
sum is a Riemann sum whose limit on m is 

lim = z\q)dq = ^-^(^(0) - zC-^)). (8) 

m-*oo zzsmU/Z Jq /smf/z ' 



We now turn to the r residual terms: Rm '■= X]j=fm/7-r+i ^^^^'^ '^^^^'^z{{2j — l)7r/m). We 
expand z around z{(2j — l)n/m) = z(2n/'j) + z'{xj) {^^^^-^^ — for some Xj £ 
(2j-i)7r 27r _ rjj_^g remainder is 

V 7 / A_/j=m/7— r+1 ' A_/j=m/7— r+1 ^7 m j ^ J ' ^ ' 

The second sum vanishes in the limit because its absolute value is smaller than '^'^"'"^^'^ sup \z!{(i) \ . 
The first sum is easily calculated and 

hm = .(^)e-^/^^^?^. (10) 
m^oo ' sm(t^/2) 

The result follows by taking the real and imaginary parts of lim(S'm + Rm) as given by ([8]) 
and ([TU]). 

(ii) The result follows from where, again, the second sum goes to zero and the factor 
z{^) = removes the r-dependent sum. □ 

Proposition 3 (Continuous case) Set 9 = ^ and denote by lim the process of taking 
the limit k,m oo while keeping 6 fixed. The thermodynamical limits o/ c'^(fci, fc2)/c| (eq. 
I^j ) with k = k2 — ki > are 

(i) (supercritical, u < Uc) 



lim e''^''d{u, g) = (1 - tanh 2z/)(l - sinh 2z/)cosh2 2z/, 

777- 

q&Qm 



independent of 9; 
(ii) (critical, u = v^) 



hm-z = 

^ sm6'/2 



(iii) (subcritical, z/ > Uc) 

lim-i Y e'''''d{u,q) 

goes to a Dirac distribution in the following sense: if f : ^ M. is a continuous 
function on the circle, then 

lim -^-YfC-^) E e-^rf(-,g) = 7/(0) 

k=l ^ ' qdQm 



with 
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2 / di^u, q) cot - dq. 
Jo 2 



We stress that only the properties stated in Lemma [T] are used in the proof and that it is 
the behavior (ii) ai q = that decides between the three regimes. 

Proof: (critical) Note that hm — i ^^^g^ e*'''^(i(z/c, g) = \im2J2q(zQ+d{i'c,q)si^kq. By 
Lemma [H d{iyc,0^) = \/2 — 1 and d{iyc,7i) = and d(uc,q) is analytic in q for q G [0,7r]. 
Therefore the limit follows from Lemma M (H)- lim —i XlgeQ™ e^'^''d{uc, q) = ^g/\ ■ 

(supercritical) To prove the supercritical case, we first replace the function d{v,q) in the 
limit by -. Then — — ^n^l. = 1 V'^'/f sm0{j 2) ^ This sum goes to ^ when m — 00 for 

any value of ^ G (0,27r). (See pLDj, App. A, Table II.) Therefore lim ^^^^^ e*«^i = \. 
Let us write q) as (i(z/, g) = ^ — + g{q) where a is the residue of diy, q) at q = 0: 
a = 2(1 — tanh 2z/)(l — sinh 2z/) cosh^ 2z/. Then the new function g is analytic on [0, tt] with 
g{0) = gi^n) = 0. This function will therefore not contribute to the limit by Lemma O (Note 
that, in the present case, there is even a supplementary factor ^ in the limit.) We thus 
obtain 

lim — - e''^''d{u,q) = - - lim — - qe''^''. 

q&Qm q&Qm 

Because | X^^^i sin 0{j — |) | < 1/| sin(6'/2) | and Yl^Ii J sin ^'(j — i) grows at most as m times 
a constant (depending on ^) as m oo, the remaining limit on the right-hand side vanishes. 
The supercritical case follows. 

(suhcritical) The limit under consideration now is limm^oo Xlfc^/ / (^) XlgeQm ^'''^^^{i', q) 
for a continuous function / on the circle T^. The case of a constant function / is simple: 

„ . m— 1 ^ . m— 1 

lim -^^/(O) J] e'^'=rf(z/,g) = /(0) lim -— ^ rf(z/, g) e^^^^ 

k=l qeQm q&Qm k=l 

= /(0) lim — V diu,q)cotl 

q&Qm 

and, because d{v, q) has a simple zero at g = 0, the summand is bounded and the limit is 
2/(0)/;rf(z/,g) cotfrfg. 

We shall now concentrate on the continuous function (7 : ^ M given by g{9) = 
f{9) - /(O) that vanishes at 6' = 0: 



^ . m-l 

lim 



Ef(^) E»"''*(''.* (11) 



m-+oo m ^-^ \ m I 

k=l ^ ^ 



q&Qn 



First note that the inner sum is an approximation of the Fourier coefficients of d{i>, g) 

1 „ X 1 r , , constant d'^d . 

- e-^rf(., q) = - e-^c^(., q)dq + -^^g^i-, 

q&Q •'^'^ ^ 
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for some g' G [— 7r,7r]. With the factor the correction terms disappear upon taking the 
hmit. Therefore the hmit ([II]) is hm„^oo -^ni Y1=q 9 (^) d-k with 4 = ^ e-*^*^rf(z/, q)dq. 
Because d{v, q) is real-analytic, its Fourier coefficients decrease exponentially, i.e. there exist 
Ci and C2 > such that 141 < Cie"'^^!^!, for /c G Z. (See, for example, V.16 of [I3-) Let 
e > and M = max^gTi |5'(6')|. Then there exists K G N such that \Y1^=k 9 (^) d-k\ < 
CiM^3^§^ < I for m > Let D = max{|(i_i|, |(i_2|, . . . , The remaining terms 

are bounded by Yl^=i 9 (^) d-k < DKmaxg^^Q 2^k-^ lfl'(^)l- Because g is continuous and 



g{0) = 0, there must be a > such that, ii m > N, then DK supg^^^ a^^xj \g{6)\ 
Therefore the limit (fTTl) vanishes and 

hm V / V e*'^'=rf(z/, q) = 7/(0) with 7 = 2/ rf(z., g) cot fdg. 

fc=i ^ ^ qeQ„ 



< I- 



□ 



2.3 The discrete and continuous cases for n < ^ 

We now turn to the general case when the number of spinflips is any number between 
and m. The calculation for n = 2 generalizes trivially to the general case and we give more 
details for this case. We drop the superindex "s" on c'^{ki,k2, ■ ■ ■ , /c2n) as it does not play 
any role here. 

To mimick the argument of the previous section, we shall write down the component, in 
(see equation ([2])), of the vector v 

1 2 3 ... ki-1 ki ki + l ... k2-l k2 ^2 + 1 ... m 
v = {- --...- + + ... + -- ... - ). 

Because the operators ^^'s (and r/^'s) are made of the operators pk and Hk that flip all spins 
to the left of the site k and may or may not change the spin k, the only vectors that may 
contribute to {^qUj)^ are the seven following (families of) configurations: 





1 




ki 




k2 


m 


I 






+ + + + + . 


. ++ 






II 


+ + + - - . 




+ + + + + . 


. ++ 






III 


+ + + + + . 


. . ++ 


+ + + + + . 


. ++ 






IV 


+ + + + + . 


. . ++ 


+ + . 


. ++ 






V 


+ + + + + . 


. . ++ 










VI 


+ + + + + . 


. . ++ 






+ + + . 




VII 


+ + + + + . 


. . ++ 






+ + + + + . 


. ++ 



The configurations I, III, V and VII have 2 spinfiips, the others have 4. The action of rjg on 
these terms in to is given by ^ie~^'^^^m~^^'^ times 
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I -e-%{h,k2) 



ki-l 



II J2 {e~''^^^'^^ - e-"'")c{l, ki, k2, m+1) 

1=2 

III (e-^9(fci-i) ^ e-i9fci)c(A;2, m + 1) 

IV J] (_e-^^('-i) + e-^«')c(fci, /, ^2, m+1) 
/=fci+i 

V (-e-^^e^^-i) _ e-'i''^)c{ki, m+1) 

m 

VI J] (e-'9(^-i) _ e-^«')c(A;i, A;2, Z, m + 1) 

l=k2+l 

VII e-^'''"c(A;i,A;2) 



Similar expressions can be obtained for the action of 77!^. We shall write the component of 
V in the equation ^gUJ — {rjg cos (f)q + rj^_g sin 0q)a; = in the form 



II + IV + VI = -( (I + VII ) + III + V). 



4 spinflips 2 spinflips 

The left-hand side is 

fcl-l fc2-l 



(+E - E + J] )e-^«'sin(0, + g/2)c(7r(/,^i,^2,m+l)) 



^=2^ /=A;i+l l=k2+l 

n IV VI 

where the symbol tt appearing in c{7i{l, ki, k2,m + 1)) is the permutation that orders the 
integers /, ki, k2, m+1. One can rewrite compactly this expression by setting c{ki, k2, k^, k^) = 
whenever two of the arguments coincide and by denoting the parity of the permutation tt 
by (-l)'W: 

m 

•gm/4gV2 ^(_iyWe-i9« sin(0^ + q/2)c{7r{l, k^, k2, m+1)). (12) 
1=2 

The right-hand side is 

e-/ieW2 + ^/2) (e-*«c(A:i, A;2) -e^^fcs^m+l) +e-^^^^c(fci, m + l)j (13) 

( I + VII ) III V 

The discrete inverse Fourier transform gives the desired expression 

(-l)'Wc(7r(A;, fci,fc2,m+l)) 



m 



J] e^'''=(i(z/, q) (e-^«^ic(A;2, m+1) - e-''^'''c{ki, m + 1) + e-^''(™+^)c(A;i, A;2)) 



The generalization to 2n spinflips is immediate. It will be obtained by examining the 
component of a (2n — 2)-spinflip configuration in {CgU)). The contribution of (2n)-spinflip 
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configurations will be of the same form as the left-hand side in (fT2|) . (Note that the above 
calculation depends only on the positions ki and k2 and on the number of stretches of 
constant signs that are flipped by the operators pk and tt/j.) Similarly the (2n — 2)-spinflip 
configurations will lead to (2n — 1) terms with alternating sign as in f|T3|) . We therefore get 
the following recursive formula. 

Proposition 4 (Recursive form) Let the number 2n ofspinflips be such that 1 <n < ^ 
and let 1 < ki < k2 < ■ ■ ■ < k2n < m be their positions. Then 

2n 

C{ki, /C2, . . . , k2n) = —J2Yl Qy''^'''~'''\-'^yc{kl, /C2, • • • , /Cj-l, kj, /Cj+i, . . . , k2n) (M) 

q&Qm i=2 

where c{ki, /c2, • • • , /^j-i, kj, /Cj+i, . . . , /c2n) = c(A;2, /cs, • • • , kj-i, /Cj+i, . . . , k2n), if n > 2 and C| 
when n = 1. 

We restrict the discussion of the thermodynamical limit to the critical regime. 

Proposition 5 (Continuous case at criticality) Set 9i = 2nki/m, 1 < i < 2n with 
n>2 and denote by lim the process of taking the limit ki,k2, . . . , k2n, m — >• oo while keeping 
the Oi 's fixed. The thermodynamical limit of c{ki, k2, . . . , k2n) is 

p{9i, 02,..., 92n) = lim —c{ki, k2, . . . , k2n) = ^^—^-TW^TT^Pi^i' O2,..., Oj, 9j+i, 62^) 

Ct sm(%i/2) 

where Oij = 6i - 6j and p(^i, 6*2, ... , 6j+i, . . . , 6'2„) = p{62, 6*3, ... , 6j^i, 6j+i, . . . , 6'2„) 

if n > 2 and = 1 when n = 1. 

Proof: This follows from Proposition [3] for u = and the elementary fact that the limit 
of a product is the product of the limits as long as these limits exist. □ 
As an example we give the continuous case for n = 2: 

P{eu 62, e„ e,) = {V2- if ( \ - \ + , ^ , ) . (is) 

\ sm 2^12 sm 2(734 sm gt/is sm 2(724 sm 2 sm 2(723 / 

This particular case was used in [2] to obtain the probability that the interface, separating 
the constant-spin clusters meeting at 6'i, ends at 62. 



3 The distribution of the number of spinflips 

Is there a "typical" number of spinflips at the boundary of a long cylinder? Or more precisely, 
what is the probability distribution of the random variable 

Y ^!L= # spinflips /2 ^ 1 
^ m ^ sites at the boundary ' 2 

as m — ^> 00? De Coninck [6] proved that the rescaled variable is Gaussian in the limit. 
We use the results of the previous section to recover this and to determine explicitly the 
mean and variance. 
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Proposition 6 The random variable = n/m on the set of configurations at the extrem- 
ity of a half-infinite cylinder behaves at criticality as 

lim Pr„ <x]=^ r e-^'/^dx 

with 11= i-^(V2 + l) anda"^ = ^(7 + 5^2) - |(3 + 2v^). 

In particular the above proposition asserts that, at the boundary of a long cylinder 
covered by a square lattice, there are on average 2m /i 2m x 0.115766 . . . spinflips, that is 
there is one at almost every 4 sites! 

Proposition [5] gives the limit of the discrete probability Prm{ki < k2 < ■ ■ ■ < k2n\n) 
for the square lattice. Unfortunately the weight included in the limit process (proportional 
to m") rules out comparing relative probabilities of all configurations with a given number 
of spinflips. One has therefore to go back to the recursive form of Proposition H] for the 
probability c*(fci, ^2, . . . , fc2n) for finite n and m. One may restrict the comparison to con- 
figurations whose spin at position 1 is + and again drop the superindex s. We shall use the 
notation d{q) for d^Ucq). We are interested in computing 

„,liSoo PT^m{Y^ = n/m) (16) 

n/m=K fixed 

where Pr^(F„ = n/m) = Ei<fci<fc2<-<fc2„<m ^(^i, ^2, • • • , /c2n)- The recursive formula ([HD 
can be used to give an explicit expression of c(fci, k2, ■ ■ ■ , k2n) as a sum over certain permu- 
tations: 

Pr™(F^ = n/m) = ^ ^ JJ d(g,-) cot (17) 

distinct g's 

We first prove Proposition [6] using (fT7j) . 

Proof: We define /(g) := d{q) cot(g) for short. It is easily checked that the above 
expression can be rewritten in the form of the probability of n successes for m/2 independent 
Bernoulh trials, each with probability pra,i ■= f{qi)/{^ + fiQi)) so 1 - pm,i = 1/(1 + /(?»))■ 
As a result the normalization is ^ = Hi^'i^l-'- + fild)- Let us write e^.i for a Bernouilli 

variable with probability pm,i- Ym is then simply ^ YIT^'^ ^rn,i- Note that the variables e^.i 
form a triangular array. (See, for example, j4|.) Moreover they are independent. Therefore 
the central limit theorem can be applied if the Lindeberg condition is verified. Since the 
variables e^.i are uniformly bounded, the condition reduces to verify that the limit of the 
variance of ^/mYm exists. But we have 

m^oo m 27r Jo (1 + f{q)) 

In addition, the mean of Ym converges to 

m/2 



1 1 r /(g) , 

m^oo m 27r Jo 1 + J (g) 
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The integrands in the expression of /i and a can be put into simple forms: 



Y^^J(^ = "^(^ + ^) (l - 2^2 + cosg + v/(l -cosg)(3-cosg) 
^ (^(-2-2^2) + (4 + 2V2)cosg-2cos^g 



+ 87r(3-2y2) 

+ (2\/2 — 2 cos q) — cos g) (3 — cos g) j . 

Integration can then be done and give /i = | — 2^("\/2 + 1) ~ 0.115766. . . and = 
(^(7 + 5^2) - 1(3 + 2v^)) ^ 0.0538198 .... □ 

We turn now to the proof of (fT7|) . Let S'2n denote the permutation groups of 2n elements. 
We shall call pfaffian a permutation vr G S2n of the integers {1,2,..., 2n} that satisfies 

7r(l)<7r(2), 7r(3)<7r(4), 7r(2n - 1) < 7r(2n) 

7r(l) < 7r(3) < ... < 7r(2n-l). (18) 

(These are the conditions on the indices of terms appearing in the pfaffian of an anti- 
symmetric 2n X 2n matrix. See, for example, [TOl E] or [H].) The set of all these per- 
mutations C S2n will be denoted Pf2n- Two well-known facts are useful here. First suppose 
that 7r(l) 7^ 1. Then there exists i > 1 such that 7r(i) = 1. If z is even, then 7r(i — 1) < 7r(i) 
and one of the pfaffian inequalities is surely false. If i is odd, then it is 7r(z — 2) < 7r(i) 
that is false. Therefore 7r(l) = 1 for vr G Pf2n- Second we obtain the cardinality of Pf2n- 
We have just seen that 7r(l) = 1. All values 2 < i < 2n are possible choices for vr(2). 
There are (2n — 1) of them. By an argument similar to the one leading to 7r(l) = 1, one 
shows that 7r(3) must be the smallest integer left in {l,2,...,2n} after deletion of 7r(l) 
and vr(2). There are then {2n — 3) choices for vr(4). Repeating the argument one gets 
|Pf2n| = (2n-l)!!(<|^2n| = (2n)!). 

Note that the previous argument for the cardinality of Pf2n follows the process by which 
the recursive expression f|T^ constructs the general term of c(/ci, /c2, . . . , k2n)- For example, 
c(/ci, ^2, ki) is given by 



— \ d{qi)d{q2){e^'^^^^^~^'^^e^'^^''^^~^^^ _ g«<7i(fci-fc3)g«g2(A:2-fc4) _j_ g«<?i(fci-fc4)g«g2{fc2-fc3)'\ _ 

gi,g2 6Qm 



The permutations (vr(l), 7r(2), 7r(3), 7r(4)) appearing here are precisely the three pfaffian per- 
mutations of Pf4. To obtain a similar form for a general n, let us denote the momentum 
introduced at the j-th use of the recursion formula by qj and the two indices on the posi- 
tions fc's by 7r(2j — 1) and vr(2j). At the first use the formula ffT^ forces 7r(l) to be 1 and 
7r(2) to be any number between 2 and 2n. This first step gives rise to the {2n — 1) terms 
(_l)ieigi(fci-fc.)_ At the second use, 7r(3) will be the smallest of the remaining integers in 
{1, 2, . . . , 2n} and 7r(4) will take any of the {2n — 3) remaining values. It is clear then that 
the indices 7r(2j — 1) and 7r(2j) appearing in ]^-|^^^.^_^ e*^^'''^''(2j-i)-^ir(2j)) qj-q those obtained by 
a pfaffian permutation of {1, 2, ... , 2n} and that all such permutations occur precisely once. 
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Because {i — 2) is the number of neighbor transpositions necessary to go from (1, 2, ... , 2n) to 
(1, i, 2, 3, . . . , i — 1, 2 + 1, . . . , 2n), the products of the factors (—1)* of ([H]) is simply (— l)'^'') 
where /(vr) is the parity of vr. One can therefore rewrite c{ki, k2, ■ ■ ■ , k2n) as follows. 

Proposition 5 (Combinatorial form) Let the number 2n of spin flips be such that 1 < 
< Y and let 1 < ki < k2 < ■ ■ ■ < k2n < m be their positions. Then 

c{k,, k2,..., k2n) = CT E (-1)'^"^ n rf(g,)e^*('-(^-^)-'=-(^^)). (19) 

The probability PTm{Ym = n/m) is then 

PTmiYra = n/ui) 

qi,q2,---,qneQm l<£<n I<fci<fc2<---<fc2n<m 7r6Pf2n l<i<" 

(20) 

Those familiar with classical works on the Ising model will not be surprised to see a pfaffian 
sum appearing here. (See, for example, |13J.) The continuous limit of (fT9|) can be deduced 
easily using the result of the previous paragraph. Using Prop. M(ii), one finds that the limit 
of the probability distribution conditionned to the number of spinflips is up to a constant 

Y n Mlio.m-(^^(2,-i))). (21) 

7rePf2n l<j<n 

This expression already exists in the literature, though in a slightly different form. Burkhardt 
and Guim [5J computed the correlation function {4>i{zi, Zi) . . . (pai^a, ^a))ci,.--,C6 fields </)j, 1 < 
i < a, in the upper-half plane when piecewise constant boundary conditions are applied along 
the real axis, with spinflips at (j, 1 < j < b. If the 0j's are taken to be the identity, one can 
argue that is nothing but the density probability that spinflips occur at ^i, . . . , (i,. 

With this identification and after a conformal map of the upper-half plane onto the cylinder 
(the point at infinity is deleted), their expression (16a) is the above expression (12T!) . Theirs 
is obtained from conformal field theory, ours from the original definition of the Ising model. 

It is interesting to remark that only the outer sum on qi,q2, ... ,qn G Qm contains the 
information about the temperature, through the function d. The rest of the expression 
is completely combinatorial in nature and rests only upon the introduction of the anti- 
commuting (fermionic) operators C,q of Section [2] in Onsager's solution. Even though they 
were introduced for the square lattice, we will see in Section H] that similar operators exist 
for the triangular and hexagonal lattices. It is therefore likely that these operators C,q may be 
introduced for a large class of two-dimensional lattices and that their commutation relations 
are independent of the lattice. The function d, on the other hand, is likely to depend on the 
lattice. It it therefore natural to introduce the function 

I<fci,fc2,.--,fc2n<ni 7rePf2„ l<j<n 
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We were not able to find any tractable form for this expression. However note that this 
expression is to appear within the sum "^^^ Yie d{Qe) where the dependency on all g's is sym- 
metric. If one decomposes Nn into its S^-symmetric component, only the fully symmetric one 
will contribute to the sum over the g's. Moreover, because the function d{q) is an odd function 
of q, it is sufficient to consider the component in Nn that is odd under the exchange of any 
of the g's. To be more specific let us introduce the following linear operators. Let P{q) be a 
polynomial in e^'^ and e~*^. The linear operator Oddg acts as (Oddg P){q) = \{P{q) — P{—q)). 
If q2, ■ ■ ■ , q-n) is a polynomial in e^^''^, e^"'^, . . . , e^**^", the linear operator Sym^^ ^j^,...,^^ 
is defined as (Sym^^ P)(gi, g2, • • • , g™) = ^ E^e5„ ?^(2), • • • , ?^_(n)) where Sn is 

the permutation group of n elements. In terms of these, the desired probability is 

Vim{Ym = n/m) = cA — \ ^ A^n(gi, 52, • • • , gn) JJ d{q() (22) 

qi,g2,---,qn&Qm i<l<n 

where 

Nn{quq2, ■■■An) 

vl<j<n / I<k\<k2<---<k2n<iri 7r6Pf2n l<i<" 

(23) 

This function Nn is defined for any positive integer n < ^ • It is a polynomial in the e^*''^ . It 
is similar in nature to the kernels of Dirichlet and Fejer arising in elementary Fourier analysis. 
Here Nn is for a discrete Fourier transform of multivariate functions whose dependency on 
their variables (momenta) is symmetric and odd. It is tempting to call it the Ising kernel. It 
turns out that this quantity has a very simple expression. It is shown in the Appendix that 
Nn{qi, q2, ■ ■ ■ , qn) vanishes whenever there exist i ^ j such that = g| and is otherwise 

^ l<j<n 

With this, the summand in VimiYrn) becomes, up to a constant, ni<j<n'^(gi) h^j when 
non-zero. This function is symmetric and even in all g's. The sum can be restricted to g's 
in = n M"*" by multiplying by 2". And because Nn vanishes whenever some of the g's 
coincide up to sign, the probability corresponds to equation (fT7|) . 

After long calculations, an agreement with explicit simulations may be sought for reas- 
surance. The Figure [2] was drawn for that purpose. The cases m = 30, 48, 60, 99, 157, 200, 
397 and 800 are plotted. They are the results of Monte-Carlo simulations except for the two 
smallest lattices that represent exact calculations. In the computer experiments the length 
of the cylinder was twice as long as the circumference and the smallest sample (5 x 10^) 
was for the m = 800 cylinder. As an example of measurements of a, this latter case gave 
a = 0.05372. The reader will notice both even and odd m's here even though the proposition 
was proved for m even only. 

The last result of this section is not a consequence of Proposition [6] as such but it does 
follow from the technique used in its proof. 
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Figure 2: Probability distributions of the variables Ym for several values of m. 

Corollary 7 The probability c-|-(m) of having only spins + at the boundary of a half-infinite 
cylinder (at the critical temperature) behaves as 

ct(m) ^ (2{V2-l)e^^/^) " (25) 



V2 

form large enough. G is Catalan's constant: G = Yl'k'=oi^^)^ ~^ 

Proof: The probability is half the probability of having no spinflip and, by the preceding 
proof, C|(m) = |Pr(n = 0) = (2 ni<i<A/(l + fM,i)Y^ with M = f and /m,* = d{qi) cot \qi, 
Qi = {2i — 1)71 /m E Q^. A direct calculation leads to 



^ _^ ^ _ 1 2^2 - 1 - cos g + ^(1 - cosg)(3 - cosg) 
sin f 2^/T^^cosq + v/2v/3~^^cosg 

Using standard tables (see for example paragraph 6.1.2 of [11}), one finds n/=i sin y 
2(m-i)/2 - Therefore 

JJ (1 + fM,i) = exp ^ log(singi/2 + d{qi) cosgi/2) - log(sin(gi/2)) 

l<i<M l<i<M 

= 2(™-i)/2exp log ^(9) 
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where 

cos q + + cos q) (3 + cos q) 
2^1 + cosg + \f2^'i + cos q 

This function /i is the remaining factor in (1 + /) evaluated at vr — g. Note that the set 
is stable under the operation g — > vr — g. Moreover h has the following useful properties: 
it is analytic on (— 7r,7r), even and takes the value 1 at g = 0. One can write log/i(g) = 
Ei>2,even ^ud therefore EgeQ+ log/i(g) = E»>2,even YJ^=\ 9j- The leading terms (in m) 
in the inner sum are g] = vr* (^2(7TT) ~ "I" ^ (m^)) • '^^^ term can be resummed 

as follows 

i>2,even ^ ^ i>2,even'^^ 

and the second as 



1 nH 71 d ^ , , , 

12 m 12m dg 

i>2,even 



TT 



24m 



The integral appearing in the first term is again somewhat difficult. Rewriting the integrand 
as follows and using Mathematica, we were able to obtain 

^ \ogh{q)dq = log (^(1 - -^)(v/l -cosg+ v/3-cosg)^ rfg = tt log(V2 - 1) + 26" 

where G is Catalan's constant. Gathering the various terms, we get (l25l) . □ 

4 Behavior at the boundary for the triangular and 
hexagonal lattices 

In this section we extend the techniques used for the square lattice to compute the quantity 
<f{ki,k2) for the hexagonal and the triangular cases. These lattices are characterized by a 
connectivity of 6 and 3 respectively except at the boundary, where more than one choices 
can be made. Our choice is drawn in figure [31 The square-shaped representation will be 
helpful for the computations of the transfer matrix that follows. 

The main result of the section asserts that the thermodynamic limit of c^{ki,k2) at 
critical temperature behaves as in the square lattice case. This yields some evidence of the 
universality of the critical behavior at the boundary. 

Proposition 8 Let m be even, k be such that 1 < k < 2m and set 6 = Let p,n and r 
as in LemmalB with 7 = 4. Denote by lim the process of taking the limit A;, m — 00 while 
keeping both 9 and r fixed. Then at criticality for k = k2 — ki 

C*(A;i,/C2) C*lattice 

lim = — — 

C| sin 9/2 

where Ciattice is either Ctn = 2(2 — v^) or Chcx = 4/3 if k is even and 2(-\/3 — l)/3 if k is 
odd. 



Even though we focused on the case of two spinflips, it is possible to generahze inductively 
our result to c^{ki, kn) in the same way as in Section 2. The result would be similar to 
the square lattice with the appropriate constant Ciattice- 

The proposition implies that, for the hexagonal lattice, various limits can be obtained 
for the same 6. To understand this, one should look at Figure [31 For the triangular lattice, 
a spinflip on an even site is actually equivalent (up to a reflection) to a spinflip on an odd 
site. This is not true for the hexagonal lattice. The case of a spinflip on an odd site is 
different from a spinflip on an even site: in the latter the sites are not nearest-neighbors 
whereas in the first, they are. This asymmetry persists in the continuum limit. A parallel 
with quantum field theory is useful. It is impossible physically to measure this probability 
precisely for a given 6. A measurement that is possible is to obtain the probability that the 
two spinflips occur within a distance 9 G [6'min, 6'max]- To calculate this probabihty (or this 
"smeared correlation function") from the above result, one would have to use the average of 
C^l^ and C^^^ as the sites with k even and those with k odd have the same density. (Of 
course this above limiting distribution is not normalizable and a cut-off becomes necessary.) 

We start by reviewing the diagonalization of the transfer matrix in Section 14.11 following 
the work of Houtappel [14] . Then we construct in Section 14.21 an equation similar to 
for the eigenvector with the largest eigenvalue. We get, in a similar way as in Section [21 an 
expression for c*(fci,/c2) which involves functions that play the same role as d{i>,q) in the 
square case. These functions are investigated only at critical temperature. As before the 
critical behavior is related to the jump of these functions at g = 0. 



4.1 Transfer matrix for the triangular and hexagonal lattices 

The element of the transfer matrices for our particular choice of lattices and boundaries for 
a transfer from row a to row a' with 2m sites (m even), i.e. the difference of Gibbs weights 
between these two rows, is 

(f 2m 2m m 

I i=i i=i i=i 

(( 2m m 
I j=l i=l 

where o"2m+i = o"i and the sum over is the sum over all possible configurations of the 
intermediate row. It is convenient to write the previous expressions in terms of matrices 
that contain interactions along columns only or along rows only to follow thereafter the 
calculation for the square lattice. We explain here the triangular case. The hexagonal case 
is done exactly the same way and only major steps are given. 

We decompose the transfer matrix from o to o' into the product of several transfer 
matrices. Figure [H depicts this process. Dotted bonds represent the identification of two 
sites (in other words a weight of 1 for the initial state and of for the other state). Plain lines 
are the usual ferromagnetic bonds. With this particular choice of decomposition interactions 
along rows and along columns lie in different matrices. One gets the following representation 
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Figure 3: Triangular and hexagonal lattices with our particular choice of boundaries. 
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Figure 4: Decomposition of the transfer matrix into Vi, V2, V3, V4 for the triangular case. 



of the transfer matrix by summing over all possible intermediate configurations 



(26) 



alalia'" 



where the Vi^s are defined as follows 



Vi(cr', 0") = exp u 



Yl ""^""j ) n ^3(fT', a) = exp I z/ ^ a,a'^ J JJ (5,^,,. 

j odd / j even \ j even / j odd 



(27) 



Second we rewrite these matrices in terms of the operators r^, k = 1, 2m and z = 1, 2, 3, 
and of the operators 7ik, pk and P (all defined in Section [2] and in [18]) 



V^i = (2 sinh 2//)^/^ exp j iv* ^ Tr^p^- j = (2 sinh 2//)"/^ exp [ iv* ^ Tr^pj j 

\ j odd / \ i even / 



2m- 1 



(2J 



V2 = V4 = exp(zi/7rip2m-P) exp -iu ^ vr^+ipj 
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where sinh z/ sinh z/* = 1. As in the square case, the parity operator P must commute with 
each Vi because the Gibbs weights are invariant under a —a. It is therefore possible to 
diagonahze simultaneously T and P. Hence we write T = |(1 + P)T^ + ^(1 — P)T~' where 
acts in a non-trivial way on the even/odd subspace and has eigenvalue zero for odd/even 
eigenvectors. Taking this into account in the previous expression of V2 and replacing Hk and 
Pk by their expression in terms of a^, we get 



j even 



Vi = (2 sinh 2//)"^/^ j _2iy* ^ a]aj - - j V3 = {2 sinh 2z/)'"/2 exp [ -2 

\ j odd / V 

= = exp ^ + ^i) 

(29) 

with T^Ij = V^^V^V^Vi and similarly for T^^.;. The condition of parity imposes a2m+i = T^i 
for V^. 

We claim that the eigenvector uj with highest eigenvalue lies in the even subspace. Indeed 
the operator P flips a configuration a into —a and it acts as — / in the odd sector. If u 
lies in the odd sector, then Pu = —u and its coefficients in the basis of configurations must 
satisfy: Co-icr2...o-2m = ~c_a-i-o-2--— 0-2^- It follows that some coefficients are negative or they are 
all zero, contradicting Frobenius theorem. We shall thus be interested in diagonalizing T"*". 

The particular form of Vi for i = 2,4 leads us to define two Fourier transforms: 

^ I ^ E,eQ„ e'^'V,, for k odd 
l$E,eQ„e^''^^.> for A; even 
r(2/-l)7r 

where Qm = i : —m 2 + 1 < / < m 2 

y 2m 

Again, m is assumed even. The particular definition of Qm allows us to invert the previous 
relations to get 

^'^ = E = E (31) 



m ^ — ' ^Jm 

k odd fc even 



The anti-commutation relations of the induce relations for the Fourier transforms 

{riq.v\} = {K^l} = ^^ {^.,^J'} = {^.,^J'} = ifg^g'. 

These relations decouple the various sectors of momentum g, — 7r/2 < q < 7r/2. As we 
shall now see the ^'s mix the sector of a given momentum q with that of momentum —q, 
but with no other. Therefore it is possible to express T in terms of a tensor product of 
operators acting on specific sectors of momenta q and —q. A direct substitution of ( 130|) in 
( |29|) gives the precise form of these operators. We get T^|; = (2 sinh 2//)™/^ no<Q<7r/2 K fo^ 
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V, = V^4T,^3,,V^2t,V^i,a where 

VS,, = exp (-2z/* (r/Jr/, + r,l^r^_^ - l) ) V;,, = exp (-2z/* (5% + - l) ) 

^2"!, = Vt,, = exp (z^e-*^ + r]-,^, + ^Jr^l^ + 5^,ri^) 

+ ('^l«^-<? + V% + (^J^g + + 

(33) 

where ^ir holds for a term identical to the preceding one but with q — > —q. A similar 
calculation of the transfer matrix for the hexagonal case yields 

W^i,, = W^3,, = exp {-2u* {{r,% + ^t^^ - 1) + ^)) 

= exp (^ue-"' [-^ivl^lg + V-A) + + + (34) 

M/4,q = exp (^z/e-^'' (-^(5jr/lg + d.^r^,) + (^Jr^, + 77I/-,)) + 

for = (2sinh2i.)™no<,<./2 for W, = W,^,W,,,Wl^W^,,. 

Recall that each Vq and Wq appearing in the tensor product are operators acting on a 
space labeled by a pair (g, —q) with q G Qm- Each space (there are m/2 in total) is of 
dimension 2^. A basis for a given space labeled by q is constructed by applying the fermionic 
operators r]±q, 5±q and their conjugates on the vacuum denoted |0). (It should be stressed 
that |0) is not an eigenvector of Vq. See |22].) One can then recognize five subspaces of a 
space q that are invariant under the action of Vq and Wq. They are spanned by: 

{^Mio)} 

{r;t|o),5t|o),5yt^t|o),5t^i^^t|o)}, 

m.5lqr^l\Q),5lr,lq\Q),5lq5l\Q),r^^^^^^^^ (35) 
{'^WJO)}. 

We will refer to these subspaces respectively as subspaces of momentum 2g, g, 0, —q and 
— 2g. Because the highest eigenvalue of the transfer matrix T is not degenerate, none of the 
highest eigenvalue of the V^'s can be either. Hence, because the actions of Vq on the subspaces 
of momenta q and —q are identical (similarly for 2q and — 2g), uj must lie in the subspace 
of momentum 0. To get cj, it thus remains to diagonalize the block acting on the space of 
momentum 0, which is of dimension 6. This can be done by writing the representation of 
Vq in this subspace in the basis fl35l) . Symbolic manipulation programs can be used to take 
exponentials and to perform the diagonalization. We will not write these expressions here 
but they will be used in the next section to compute d'iki, k2). 

4.2 Computation of c^(A;i, A;2) 

The first step in computing c^{ki, /C2) is to construct two operators that annihilate u, instead 
of one for the square lattice. It then remains to apply them to u in the basis of the spinflips 
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to get c^{ki, ^2)- We will need three important facts. The results apply to both the hexagonal 
and triangular cases. Vg stands here for either Vg or Wg. 



Lemma 9 For all q > in Qm, the operator Vg satisfies the relation <p'^Vj'(j)gVg = 1 where 
(j)q = (5j + 6g) + (r^J + rig) (77!^ + rj^g) and " " denotes the transposed operator. 

Proof: By definition of 5g and rjg in terms of Pauli matrices, one gets 0^ = (pg and 
0^ = 1. Hence it suffices to show that = 1, = 1, 2, 3, 4. This is directly verified 

by expressing every operator in terms of 5g and rjg and by using their anti-commutation 
relations ([32]) • □ 

Proposition 10 Let uoi and uj2 he two eigenvectors of Vg with eigenvalues Ai and A2 re- 
spectively. Then either Ai = I/A2 or {011,002)^^ = 0, where ( ■ , ■ is the bilinear form 
induced by (j)q: {u,v)^^ = u'^cpgV. 

Proof: By direct calculations, using Lemma [9l □ 

Corollary 11 Ifvg is an eigenvector ofVg with eigenvalue A 7^ 1, then {vg,Vq)^^ = 0. 

We note that the subspace of momentum is an invariant subspace of (pg so that the 
corollary holds on this subspace. If Ug denotes the component of co in the subspace q of 
momentum 0, we must have by the corollary that {ug,Ug)(f,^ = 0, since the action of Vg and 
Wg on the subspace is different than the identity. In the basis of ( !35|) . this becomes for 

ViVe - V2V3 - = 0. (36) 

We stress that all the above statements hold at any temperature. 

The desired operators are now constructed. By analogy with the square case we choose 
the form 

= ar]g + b6g + cr]lg + d5^_g. (37) 

where a, 6, c, d may depend on q. We argued previously that u belongs to the subspace of 
momentum 0. Hence, from the definition of ^g, ^gUJg must lie in the subspace of momentum 
—q. The equation ^gUjg = written in the basis becomes in this 4-dimensional subspace 

/ V5 -vs -vi \ /«\ 

-f 2 -f 4 f 1 ^ - 

—Vq f4 f3 C 
\ Vq V2 -V5J \dj 

If we assume that the components Vi of u are non-zero, then the 4x4 matrix above can be 
shown to be of rank 2 by fl36|) . It is easily checked that the kernel is spanned by (— f2, f5, 0, vq) 
and (f4, — f2, Vq, 0). These vectors yield the two desired operators that we now refer to as ^g 
and ^'g respectively. 

The coefficient c^{ki,k2) can now be calculated as in section [2J We first write uo in the 
basis of the spinfiips: uj = c^a^ + qa^ + ^"^^ Ei<fei<...<fc2;<2m ^20(±; •-, k2i). 

We compute the action of ^g and ^g, < q < 7r/2, on uj in the space of momentum taking 
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and in the form (137|) . We recall that the action of these operators is well understood 
in the basis of spinflips if one writes and in terms of operator Tr^ and pk that change 
the sign of the spins from the site 1 to — 1 and from 1 to respectively. The result of the 
action must be 0. We calculate c^{l,k) and c^(fci,/c2) follows by translational invariance. 
The projections of and ^'gU on yield two equations: 



(3J 



ct {-vse -v^- iv^) = ^ c (1, k)e '"^ {-v^e''^ + - iV(^) 

k even 

+ ^ c-(l, k)e-''^^ {v, - v.e''^ - tvee''') 

k odd 
k even 
k odd 

These two equations can be solved for c~(l, k) and ^^^^^ c~(l, A;). The inverse 

Fourier transform on q yields c~(l, k). By translational invariance, we get for any ki, k2'- 

,s^^ ^^ f2^E,eQ.e^^^'^-'^)4ven(^,g), if fc^ - ^1 + 1 is evcu, 

l2^E,eQ.e^''^^-^'^^^odd(^,g), if A;2 - A;i + 1 is odd, ^ ^ 

where 

rfevenl^^, 



c^oddl^^, q) 



sinqivQ - - cosq^v^ + Vr,) + (vg - 172)' 
2 cosg(l - + 2sing(t;4 + v^) + 2i(t;2 + fs) 



sin g(t;6 - Vi) - cos g(f 4 + f 5) + (t^s - ^2) 
The components fi,...,f6 of cUg are obtained following Section [4.11 The expressions for a 
generic temperature v are rather heavy. We write down the result for v = v^.. The use of 
symbolic manipulation programs was essential to obtain the expressions in these relatively 
simple forms. 
Triangular Lattice 

{2(-3-sing+2(cos l+sin 1) v^3+sIrrg)4-2i(2cos ij-(cos |-sin 2) v^3+sIFg) -fn^ ^ it 
-Re Kven(i^c, -q)] + t Im Kven(^c, -q)] , if -f < g < 0; 

4(^cos |-sin 1) (^cos l+sin l-iv^S+smg) ■( n <^ <^ i" 

Cd(^c, q) = { ^ ' ^ ^ ^ ^ 2 ; 



(40) 



Hexagonal Lattice 



2(v^— 1) (s+sing-VS^cos 2+sin l^ys+img) -f n <^ ^ 

sinn(3+sinQ) ' II U < g < ^ ; 



Ce„(^c,g)=<^ ^ ^ sin,(3+sin,) 

-C^ev'enl^^c-g), if-f <g<0; 

I (cos f-sin 1) (-2v^V3+sin<j+2(cos |+sin |)(3+sin<j)) --f n ^ ^ 

rfhex (^^, g) = — — ^ sin,(3+.n,) ' it U < g < ^; 

.-<d(^c,-g), if-f <g<0. 
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(41) 



We gather in a lemma the properties of the above functions that are relevant for the limit 
m oo. Its proof is similar to that of Lemma [H The continuum limit of 2mc'^{ki, k2)/c^ in 
Proposition [S] follows directly from Lemmas [2] and \n\ with rO = 2iTk/'~f = nk/2. 

Lemma 12 The four functions d^Jdd' ^even; ^odd? '^even [~f'f] have the following proper- 
ties: 

(i) their restriction to [0, ^] satisfy the requirements put on z in Lemma^ 

(ii) their parity and some of their relevant values are as follows. 





values in 




parity 


(^^c,0+) 


(^c, f ) 


Jtri 
"odd 


M 




odd 


2(2 - V3) 





even 


C 


tDq Jtri 

"even 

T Jtri 
llli Ugygj^ 


odd 
even 


2(2 - V3) 
2(2 - V3) 





^hex 
"odd 
Jhex 
even 


M 
M 




odd 
odd 


4 

l(v^'-l) 






We were not able to bring the expressions for a generic temperature u into a form suitable 
for publication. Like the function d of the square lattice (eq. (jl])), they seem to have an 
interval of discontinuity in the (z/, g)-plane extending from (0, 0) to (z/c, 0). This is shown on 
Figure [5] where the real part of d^J^^^ has been plotted. Because of the singularity along the 
interval with q = and u G (0, z/^], we have to limit the range (vertical axis) to the interval 
[—0.8, 0.8]. The cut on the surface imposed by this limitation has been marked by a darker 
curve. The restrictions of Re eleven (^^^ Q') three values of u have been also highlighted. 
These three values are J^c + ^ = 0.299653 ... in the subcritical phase, the critical value 
z/c = 0.274653 . . . and J^c ~ ^ = 0.249653 ... in the supercritical phase. The subcritical 
curve is obviously smooth. Lemma [T2] has established that the cricital curve is real-analytic 
except at g = where there is a jump: Re d^^^^^(uc,0^) = —Re (ig"g„(z/c, 0~) = 2(2 — v^). 
The supercritical curve seems to have a pole at g = 0. 



Appendix: A combinatorial lemma 

The purpose of this appendix is to show the following lemma. 

Lemma 13 (Ising kernel) Let m be a positive even integer and n a positive integer such 
that 2n < m. The polynomial Nn{qi, q2, ■ ■ ■ ,qn) in e^*''^, e^"'^, . . . , e^**^" defined by 

n OddJ Sym,^,,^,„„,„ Yl E (-1)'^^^ n e^(^^<^-^'-'-<^^')* (42) 

is 

if there exists i ^ j such that qf = g|, 

TT 1 r 

II cot^gj, otherwise, 
i<i<n 

when the values of the q's are restricted to the set Qm = {(2j — l)7r/m, — y + 1 < j < y}. 
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(subcriticni) 



Figure 5: The function Re eleven (^^ ^) for the triangular lattice. 

The notation is that of Section [31 
Proof: We first establish a one-to-one correspondence between summands in 

I Y[ OddqASym^^ ^j^^ ^ (-l)'(^)p(g^(i), g^(2), . . . , g^(„)) 

\l<j<n / 7rePf2„ 

and in 

J2 i-^y^'^Piqpii),gp{2),---,qp{n))- 

Take any permutation p in S'2„. There exists a unique choice of n permutations G S2n, 
acting trivially on all indices but 2i — 1 and 2i such that ri(p(l)) < ri(p(2)), r2(p(3)) < 
T2(p(4)), Tn{p{2n — 1)) < r„(p(2n)). After the pairs have been ordered, there is 

a single permutation a G that orders pairs in increasing value of their first elements: 
cr(ri(p(l))) < (T(r2(p(3))) < ■ ■ ■ < a{Tn{p{2n — 1))). Because the choice of r^'s amounts to 
choose one of the two terms in Oddg- and the choice of a is one of the terms of the sum 
Sym^^ g2,...,g„) one-to-one correspondence stated above follows with aTiT2 . . . r„p = tt. In 
Oddq;, the term in which q is replaced by — g is multiplied by —1 and the corresponding 
term in ^^g^^ should appear with a factor Hil"-'-)'''^'''- There is no alternating sign in the 
operator Sym. However, because a permutes pairs, it is an even permutation. Therefore 
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ll.{-iy(-^) = {-iy(p^ and 



iV„(gi, g^, . . . , g„) = ^ Yl E (-1)'^'^ n e^('=''(^-^'-'''(^^')*. (43) 

Note that by the introduction of auxihary variables X2j-i = e*'^^ and X2j = e~** for 1 < j < n, 
the sum over the permutation group S2n becomes simply a determinant: 

I<fci<fe2<---<fc2n<m. 

The second step is to transform the sum over ordered 2n-tuplets {ki,k2, . . . , k2n) into a 
sum over partitions and to introduce the appropriate tools from the theory of symmetric 
functions. A partition A is a finite set A = (Ai, A2, . . . ) of positive non- increasing integers. 
(See p5j for standard notations and definitions on partitions.) Define the integers Aj, 1 < 
i < 2n hy Xi = k2n-i+i — {2n — i + 1) or kj = j + X2n+i-j- The ordering 1 < fci < ^2 < 
■ ■ ■ < ^2n < is equivalent to m — 2n > Ai > A2 > . . . A2n > 0. The partial ordering fi G X 
between partitions stands for /ij < Xi, for all i. Finally (M^) stands for the partition 

(M^) = ( M,M,...,M ,0,0,...). 

N times 

With this notation, the polynomial Nn can be rewritten as 

Ac((m-2n)2") 

One can change the order of rows (1,2,..., 2n) to {2n, 2n — 1, . . . , 1) in the determinant. The 
number of transpositions has the parity of (2r;,)(2n — l)/2 = n{2n — 1) and has therefore the 

parity of n. Then det(a;-"^^''"+'"')i<ij<2„ = (-1)" det(a;f-''*"^"~^)i<ij<2„ ni<j<2n ^i- the 
present case, the value of the x/s are such that X2j-i = l/a;2j and the last product is simply 
unity. The definition of the Schur function of N variables associated with the partition 
A (also to be found in [TBJ) is sx{xi,X2, ■ ■ ■ , xjy) = ax+s{xi,X2, ■ ■ ■ , x^) / as{xi,X2, ■ ■ ■ , x^) 
with ax+sixi, X2, ■ ■ ■ , xn) = det{xf^^^ '')i<i,j<N and as is the Vandermonde determinant 
as{xi,X2, . . .,xn) = det(xf"^)i<ij<Ar = Yli^i^j^^ixi - xj). The polynomial Nn is then 

(—1)" 

2^2; • • • 5 X2n)- (46) 

i<i<j<N Ac((m-2n)2") 

This form has the remarkable advantage over previous ones that the above sum over parti- 
tions is known to combinatorists. On page 84 of [15] one finds X]ac(a/^) ^xi^i^ X2, ■ ■ ■ , xn) = 
Dm/Do where Dm = det{xf™~' - ^)i<*j<iv and Do = det(a;^^^"* - xf-%<ij<N. The 
computation of Nn is then simply that of two determinants. 
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The third and last step is the computation of Dq and Dm- Both shares some obvious 
factors. Note that M + 2N — m + 2n and N — 2n and their forms are 



^0 



An-2 



X 



- 1 

- Xx 



xT~^ - 1 

^4n 2 — 



X 



4n-l 
2n 



- 1 



X 



4n-2 
2n 



— X2n 



X 



2n 



X 



2n-l 



X. 



.2n 



X. 



2n-l 



X 



2n 
2n 



X 



2n-l 
2n 



and 



D 



M 



^m+2n-l _ 
m+2n-2 



X 



X\ 



- Xi 



^m+2n-l _ ^ 
^l+2n-2 _ 



X. 



m+2n—\ 
2n 



- 1 



X\ 



m+2n-2 
2n 



- X2r. 



X 



2n-l 
1 



X. 



2n-l 



X2n X 



2n-l 
2n 



Each element of the j-th row of both Dm and Dq has an obvious (xj — 1) factor. (Recall 
that m > 2n and that m + 2n > 2i for 1 < i < 2n.) Both determinants must contain a 



l<j<2n\-^J 



factor Yl 

contain abactor Ui<i<j<2n^ 
Dm, it becomes 



{xj — 1). If Xi = Xj for i j, these determinants vanish. They must also 
- Xj). Finally, if one does Xj — 1/xi in the j-th column of 



(j-th column)!^. 



= 1/Xi 



-X, 



—■m—2n+k 



Xi )l<k<2n 



-m—2n+l ( ^m+2n—k 



X, 



— X' 



h<k<2n 



= -X, 



-m—2n+l . 



i-th column) 



(47) 

(48) 
(49) 



Ll<i<j<2n 



{xiXj — 1). A similar computation for Dq shows 



and Dm must contain also a factor 
that it shares also this factor. 

Let us define the leading monomial in a polynomial P G Q[a;i, 2:2, . . . ,a;2rt] as the one 
obtained through the following procedure. Among all monomials of P find those that have 
the largest exponent of Xi. Among the latter find those that have the largest exponent of 
X2- And so on till one chooses the one with the largest exponent of X2n- This monomial 
is the leading one. The leading monomial in D^ is easy to determine. From the above 
determinantal form, the largest exponent of Xi is 4n — 1. Deleting the first column and row, 
we conclude that the largest exponent of X2 among the monomials containing x\"'^^ is An — 2. 
And so on. The leading monomial of Dq is Y[i<j<2n^T~^ coefficient is unity. The 

leading monomial in the factors identified above in Dq 



\<j<2n 



.-1) n 

l<i<J<2n 



Xj)(XiXj 1) 



can be also determined. The first product contributes 1 to the exponent of Xi in the leading 
monomial and the second 2{2n — 1). This exponent is therefore (4n — 1). Among those 
monomials containing the first product contributes 1 to the exponent of X2 and the 

second (2n — 2) + (2n — 1) and the exponent of X2 in the leading monomial is 4n — 2. 
Repeating the argument, it is clear that the leading monomial in the factors above coincide 
with that of Dq and its coefficient is also 1. We have thus shown that Dq = do and verified 
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that Dm has Dq as factor (but this is a basic fact in the theory of Schur functions). Thus 
dm = Dm/ Dq is a polynomial G Q[xi, X2, . . . , X2n\- The polynomial has still a new form 
Nn{qi, q2, ■ ■ ■ , Qn) = ^"^^CLsdrn- This simple form has as immediate consequence that iV„ 
vanishes whenever Xi = xj with i ^ j, that is when there exists a pair 1 < A;, / < n of distinct 
integers such that qk = qi {i and j have the same parity) or qk = —qi {i and j have different 
parity). 

To calculate the value of A^„ when all qf are distinct, let us first remark that some of the 
factors in as and in Dq give easily the factors cot of the proposed answer. Indeed the 
quotient of the factors ni<j<n(^2j-i — X2j) in as with the factors ni<j<2n(^i ~ -'-) -^0 is 



n,<j<„(e-" - - 1) V (2i)V iism^^,,. ^J-l„ 2 

Let us denote by (i, j) pairs of odd integers with 1 < i < j < 2n. Then the residual factors in 
as, that is those not used in the above quotient, are Y\^^ ■^{xi — Xj){xi+l — Xj){xi — Xj^l){xi+l — 
Xj+l) and those in Dq are ni<i<j<2n(^* ~ ~ 1) 

^ (50) 



Ul<j<ni^2j~l - X2j){x2j-lX2j - 1) U{i,j){XiXj - l){Xi+lXj ~ l){XiXj+i - l)(Xi+iXj+i - 1) ' 

The strategy will be to first simplify the factors {x2j-i — X2j){x2j-iX2j — 1) in Dm and then 
use freely the relation X2j-i = X2^ in the remaining expression to identify the last factors of 
dSOD. 

The factors {x2j-i — X2j){x2j-iX2j — 1) occurs in any 2x2 determinant of the elements in 
the rows {2k — 1) and {2k) and in the lines i and j, i < j. If we set x = X2k-i and y = X2k, 
this is 



dij{x,y) 
A direct calculation gives 



^m+2n-i x^~^ ym+2n-i y^~^ 

^m+2n-j _ ym+2n-j _ yj~^ 



(X 



m+2n—2—i—j 

- y){xy - 1) ( {x/yf ^ {xy^ 

k=0 1=0 

m+2n—2—i~j j — t— 1 



_ ^m+2n-2-j 



fc=0 /=0 



After simplification with the factors {x — y){xy — 1) = {x2j^i — X2j){x2j-iX2j — 1) of ([50 
we can then use the identity x = y^^ in the remaining expression dij {x, y)/{x — y) {xy — 1) ■ 
ujij{x)/{x — x~^) where 

UJij{x) = (m + 2n - 1 - i - j){x^-' - X'-^) - {j - z)(x"+2n-l-i-i _ ^-{m+2n-l-i^j)y 

29 



Define the variables yj = X2j-i (and y- ^ = X2j) for 1 < j < n. The remaining factors of (150!) 
are 

JJ(a;iXj - l){xi+iXj - l){xiXj+i - l){xi+iXj+i - 1) 

= n (51) 

l<i<j<n 

and the determinant Dm has the following form in terms of the Uij-. 

LLl<j<nyyj yj ) nePhn l<j<n 

(Note that there are {2n — 1)!! permutations in Pf2„ and that the symmetrization operator 
gives rise to n! terms out of each summand. The right-hand side contains (2n — l)!!n! = 
(2n)!/2" terms. This is the right number as each of the n factors cu's is itself the sum of two 
terms.) We are now interested in the value of Dm at y/s such that none of the factors in 
the denominator is zero. We can therefore evaluate Dm (and the cj's) at yj G Qm, that is we 
can use in Dm the fact that yj^ = —1. (This is the first time since A^^ has been expressed in 
terms of the quotient Dm/Dq that this is possible.) The cu's are 



u;ij[x) 



{m + 2n-l-i- j){x^-' - x'-^) + (j - i){x^''-^-'-^ - a;-(2"-i-»-J> 



for X a m-root of —1. We can determine the leading monomial of the polynomial Dm (reduced 
by replacing all yj^ by —1) for the variables 1/1,1/2, ■ ■ ■ ,yn taken in that order, yi being the 
first. Because 1 < i < j < 2n, the first term of any Uij reaches its largest value {2n — 1) 
when i = 1 and j = 2n and the other terms have smaller exponents than this maximum. 
The leading term in iOi^2n{yi) is then my\^~^. The range of i and j for tOij{y2) is restricted to 
2<i<j<2n — l. An analysis similar to the one just done reveals that i = 2, j = 2n — l are 
the values to choose and the corresponding u has leading term my"^'^ . This process leads 
to the leading monomial in the above sum over vr G Pf2n- It is Hj Z/j""*'^"'"^''- Its coefficient is 

m". The corresponding monomial for Dm is therefore Y\ - y^*"" ''^ with the same coefficient. 
(Note that the parity of the permutation putting (1, 2n, 2, 2n — 1, . . . , n — 1, n) in increasing 
order is even and that it is Pfaffian.) 

Similarly, for the largest exponent of yi in ( 15T|) . the first parenthesis contributes (n — 1), 
the second also (n — 1) and the two last ones contribute nothing to the exponent but gives a 
factor (— iq i|;g coefficient. The exponent is then (2n — 2). Repeating the argument for 
the other variables, we find that the leading monomial is the same as for Dm and its coefficient 
is unity. We must therefore conclude that, upon evaluation of the tu's at m-roots of —1 whose 
squares are distinct, the quotient of Dm with the residual factors (|5U]) is independent of the 
qj chosen and is rrf. The polynomial Nn takes then the value ^(— ^)" IlKjXn cot |gj 
whenever the squares g| are distinct. □ 

One has to note that the last calculation that rests upon the use of yj G Qm leads to 
the wrong answer whenever yi = yj or y^ = y~^ for a given pair of distinct integers i and 
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j. (The value of Nn at these points was obtained earlier in the proof.) This wrong result is 
not surprising as both Dm and (1511) vanish at these points. One has to cancel out common 
factors before using the property yj G Qm- We have done this exercise only for m = 4, n = 2 
which is the first non-trivial case. In addition to the m" = 4^ terms computed above, one 
gets the following polynomial, after simplification of the common factors and then use of 
Vj e Qm- -8 - 4:xy + Ax^y + Axy^ - x^y^ which, for x, y G Qm, is -niJ^{5^y^i + It is 

obviously the right correction for the theorem to hold. 
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